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We show that the number of fc-matching in a given undirected graph 
^ ' G is equal to the number of perfect matching of the corresponding graph 

Gk on an even number of vertices divided by a suitable factor. If G is bi- 
partite then one can construct a bipartite Gk- For bipartite graphs this 
^ ■ result implies that the number of fc-matching has a polynomial-time ap- 

. proximation algorithm. The above results are extended to permanents 

and hafnians of corresponding matrices. 
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1 Introduction 

5r ! Let G = (V, E) be an undirected graph, (with no self-loops), on the set of 

vertices V and the set of edges E. A set of edges M C E is called a matching 
if no two distinct edges 61,62 € M have a common vertex. M is called a 
k-matching if #M = k. For G N let M.k{G) be the set of /c-matchings in 
G. {A4k{G) = for A; > [^J-) If #^ = 2n is even then an n-matching is 
called a perfect matching. (j){k, G) := j^A4.k{G) is number of /c-matchings, and 
let 0(0, G) := 1. Then $(x, G) := XlfcLo '^(^' C!)x^ is the matching polynomial 
of G. It is known that a nonconstant matching polynomial of G has only real 
negative roots [H]. 

Let G be a bipartite graph, i.e., V = Vi U V2 and E C Vi x V2. In the 
special case of a bipartite graph where n = ^Vi = #V2, it is well known that 
(j){n,G) is given as perm B{G), the permanent of the incidence matrix B{G) 
of the bipartite graph G. It was shown by Valiant that the computation of the 
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permanent of a (0, 1) matrix is #P-complete jHj. Hence, it is believed that the 
computation of the number of perfect matching in a general bipartite graph 
satisfying = cannot be polynomial. 

In a recent paper Jerrum, Sinclair and Vigoda gave a fully-polynomial ran- 
domized approximation scheme [fpras) to compute the permanent of a non- 
negative matrix [7|. (See also Barvinok [T] for computing the permanents 
within a simply exponential factor, and Friedland, Rider and Zeitouni jS] for 
concentration of permanent estimators for certain large positive matrices.) 

[7j yields the existence a fpras to compute the number of perfect matchings 
in a general bipartite graph satisfying j^Vi = #l2- The aim of this note is to 
show that there exists fpras to compute the number of fc-matchings for any 
bipartite graph G and any integer A; G [1, In particular, the generating 
matching polynomial of any bipartite graph G has a fpras. This observation 
can be used to find a fast computable approximation to the pressure function, 
as discussed in jl], for certain families of infinite graphs appearing in many 
models of statistical mechanics, like the integer lattice Z"^ . 

More generally, there exists a fpras for computing perm^ B, the sum of all 
k X k subpermanents of an m x n matrix B, for any nonnegative B. This is 
done by showing that perm^, B = i^„T^l^(^n-ky. ^ corresponding (m + n — 
k) X [m -\- n — k) matrix Bk- 

It is known that for a nonbipartite graph G on 2n vertices, the number of 
perfect matchings is given by hafyl(G'), the hafnian of the incidence matrix 
A{G) of G. The existence of a fpras for computing the number of perfect 
matching for any undirected graph G on even number of vertices is an open 
problem. (The probabilistic algorithm suggested in j7j applies to the compu- 
tation of perfect matchings in G, however it is not known if this algorithm is 
fpras.) The number of fc-matchings in a graph G is equal to h.a.fkA{G), the 
sum of the hafnians of all 2k x 2k principle submatrices of A{G). We show 
that that for any m x m matrix A there exists a (2m — 2k) x {2m — 2k) ma- 
trix Ak such that haf^ A = (2m^fc)! • Hence the computation of the number of 
/c-matching in an arbitrary G, where n = 0{k), has fpras if and only if the 
number of perfect matching in G has fpras. 



2 The equality perm^ B = 

Recall that for a square matrix A = [aij]^j=i G M*^^*^, the permanent of A is 
given as permy4 := J2aeSn '^io-(i) • • • '^na{n), where Sn is the permutation group 
on (n) := Let Qk,m denote the set of all subset of cardinality k of 

(m). Identify a G Qk,m with the subset {ai, . . . , ak} where 1 < ai < . . . < 
ak < m. Given an m x n matrix B = [hjJ^jZi ^ M"*^"- and a G Qk,m, P ^ Qi,n 
we let B[a, (3] := [baii3j]ifj=i G M'^^' to be the corresponding k x I submatrix of 
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B. Then 

peim^B := ^ peTmB[a,P]. 

"6Qk,m,/36Qk,n 

Let G = (Vi U V2, E) be a bipartite graph on two classes of vertices Vi 
and V2. For simphcity of notation we assume that E G Vi x V2. It would be 
convenient to assume that = 171,^X2 = n. So G is presented by (0,1) 
matrix B{G) E {0, l}'"^". That is B{G) = [bij]'^!l-^ and kj = 1 ^ (ij) E 
E. Let k E [1, min(m, ri)] be an integer. Then /c-matching is a choice of k 
edges in Ek := {ei, . . . , Ck] C i? such that Ek covers 2k vertices in G. That is, 
no two edges in Ek have a common vertex. It is straightforward to show that 
perm^, B{G) is the number of A;-matching in G. 

More generally, let B = [bij] E M'"^", M+ := [0,oo) be an m x ri non- 
negative matrix. We associate with B the following bipartite graph G{B) = 
{Vi{B)UV2iB),E{B)). Identify Vi{B),V2{B) with (m), (n) respectively. Then 
for i E {m),j E (n) the edge is in E{B) if and only if bij > 0. Let 

Gu, := (V^i(-B) U V2{B),E^{B)) be the weighted graph corresponding to B. 
I.e., the weight of the edge E E{B) is bij > 0. Hence B{Gw), the rep- 
resentation matrix of the weighted bipartite graph Gyj, is equal to B. Let 
M E M.kiG{B)). Then n(.ij)eAf^ji weight of the matching M in G^^,. 

In particular, perm^ B is the tota/ weight of weighted /c-matchings of Gw The 
weighted matching polynomial corresponding to B E M™^", or G^^, induced by 
S, is defined as: 

min(m,n) 

S) := perm^ B x\ B E M+^", permg S := 0. 

fc=0 

can be viewed as the grand partition function for the monomer-dimer 
model in statistical mechanics (See §3 for the case of a nonbipartite graph.) 
In particular, all roots of $(x, B) are negative. 

Theorem 2.1 Let B E and k E (min(m,n)). Let 

Bk E 5g ifif, following 2x2 block matrix 

Bk := ^ ^ '^m,rn-k ^ whcrc Ip^q is a p X q matrix whose all entries are 

equal to 1. Then 

perm Bk 

P^^"^^'^= (m-fc)!(n-A:)! - ^^''^ 

Proof. For simplicity of the exposition we assume that k < min(m, n). (In 
the case that k = min(m, n) then Bk has one of the following block structure: 
1 X 1, 1 X 2, 2 X 1.) Let G^ = (Vi(5) U ^2(5), (5)), G^,^ = (V.iBk) U 
V2{B), E^{Bk)) be the weighted graphs corresponding to B,Bk respectively. 
Note that G^ is a weighted subgraph of G^^k induced by Vi{B) = (m) C 
{m + n — k) = Vi{Bk), V2{B) = {n) G {n + m — k) = V2{Bk). Furthermore, 
each vertex in Vi{Bk)\Vi{B) is connected exactly to each vertex in V2{B), and 
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each vertex in V2{Bk)\V2{B) is connected exactly to each vertex in Vi{B). The 
weights of each of these edges is 1. These are all edges in G{Bk). A perfect 
match in G{Bk) correspond to: 

• An n — k match between the set of vertices Vi{Bk)\Vi{B) and the set of 
vertices /?' G Qn-k,n, viewed subset of 1^2(5) . 

• An m — k match between the set of vertices V2{Bk)\V2{B) and the set 
of vertices a' G Qm-k,m, viewed subset of Vi{B). 

• A k match between the set of vertices a := (m)\a' C Vi{B) and f3 : = 
{n)\(3' C V2{B). 

Fix a G Qk,m, f3 G Qk,n- Then the total weight of /c-matchings in Gw{Bk) 
using the set of vertices a C V\{Bk),P C V2{Bk) is given by perm 
The total weight of n — A; matchings using Vi{Bk)\Vi{B) and (3' C V2{Bk) 
is [n — k)\. The total weight oi m — k matchings using V2{Bk)\V2{B) and 
a' C Vi{Bk) is {m—k)\. Hence the total weight of perfect matchings in G^{Bk), 
which matches the set of vertices a C Vi{Bk) with the set [3 C V2{Bk) is given 
hy {m — k)\{n — k)\'geim.B[a, 13] . Thus perm 5^ = (m — /i;)!(n — /c)! perm^ S. □ 

We remark that the special case of Theorem 12.11 where m = n appears in 
an equivalent form in j2] . 

Proposition 2.2 The complexity of computing the number of k-matchings 
in a bipartite graph G = {Vi U V2, E), where 

min(#Vi,#V2) >k> cmax(#Vi, #V2)" and c,a G (0,1], is polynomially 
equivalent to the complexity of computing the number of perfect matching in a 
bipartite graph G' = {VI U V^, E'), where #V/ = #1^2' ■ 

Proof. Assume first that G = {ViU V2,E),m = #14, n = #^2 and 
k G [cmax(#Vi, #1^2)°, inin(m, n)] are given. Let G' = (y(UV2, E') be the bi- 
partite graph constructed in the proof of Theorem l2.1l Theorem 12 . II vields that 
the number of perfect matching in G' determines the number of A;-matching 
in G. Note that n' := = #^2' — 0{ka). So the A;-matching problem is a 
special case of the perfect matching problem. 

Assume second that G' = {V( U V2, E') is a given bipartite graph with k = 
^Vi = #^2- Let m,n>k and denote by G = {V1UV2, E'),i^Vi = m, #1/2 = n 
the graph obtained from G by adding m — k,n — k isolated vertices to V(, V^' re- 
spectively, {E' = E) . Then a perfect matching in G' is a /c- matching in G, and 
the number of perfect matching in G' is equal to the number of A;-matchings 
in G. Furthermore if A; > cmax(m,ri)" it follows that m,n = 

The results of [7j yield. 

Corollary 2.3 Let B G M!j!^" and k G (min(m,n)). Then there exists a 
fully-polynomial randomized approximation scheme to compute perm;. B. Fur- 
thermore for each x G M there exists a fully-polynomial randomized approxi- 
mation scheme to compute the matching polynomial B). 
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3 Hafnians 



Let G = (y,E) be an undirected graph on m := vertices. Identify V 

G {0, 1}™^™ be the incidence matrix of G, i.e. 



with (m). Let A{G) 

Qij = 1 if and only if (i, j) G E. Since we assume that G ia undirected and has 
no self-loops, A{G) is a symmetric (0, 1) matrix with a zero diagonal. Denote 
by Sm(T) D Sm^o('^) the set of symmetric matrices and the subset of symmetric 
matrices with zero diagonal respectively, whose nonzero entries are in the set 



T C M. Thus any A = [aij] e Sm,o(K+) induces G{A) = 
V{A) = (m) and {i,j) G E{A) if and only if Uij > 0. 



: {V{A),E{A)), where 
Such an A induces a 



weighted graph (^^(A), where the edge G E{A) has the weight aij > 0. 
Let M G Mk{G{A)) be a A;-matching in G{A). Then the weight of M in 
G^{A) is given by n(ij)6M«ii- 

Assume that m is even, i.e. m = 2n. It is well known that the number 
of perfect matchings in G is given by hafyl(G), the hafnian of A{G). More 
general, the total weight of all weighted perfect matchings of Gw{A),A G 
S2n,o(I^+) is given by haf A, the hafnian of A. 

Recall the definition of the hafnian of 2n x 2n real symmetric matrix A = 
[aij] G M^*^^^". Let i^2n be the complete graph on 2n vertices, and denote 
by Ai{K2n) the set of all perfect matches in K2n- Then a G A4{K2n) can 
be represented as a = {12-, 32), ••, (^n, jn)} with ik < jk for A; = 1, . . .. 



Denote 



=1 '^ikjk- 



Then haf A 



€M{K2n) 



Note that haf A 



2nx2n 

nxn -g j-gia^gfj the determinant 



does not depend on the diagonal entries of A. Hafnian of A is related to the 
pfaffian of the skew symmetric matrix B = [ba] G M^"^^", where = a,;,- if 
i < j, the same way the permanent of C G K 
of G. Recall pfaf 5 = ^^^^^^^^^ sgn(a)6„, 

the permutation a G S2n given by a = 



det5 = (pfaf S)2 
Let A G S 



where sgn(a^ 
2 3 4.. 
ii i2 32 ■■ 



is the signature of 
2n 



Jr. 



Furthermore 



Then 



haffcy4:= haf74[a,a], k = 1, . . . , \_—\. 



For A G Sm,o(IR+) haffc74 is the total weight of all weighted fc-matchings 
in Gw{A). Let hafo(A) := 1. Then the weighted matching polynomial of 
Gyj{A) is given by ^4) := 'Yl,k=o ^^^k A x^. It is known that a nonconstant 
^4), A G Sm,o(lR+) has only real negative roots [Oj. 

Theorem 3.1 Let A G 5m,o(K+) and k E (Lf J)- Let Ak G S2m-2k,o(K+) 
be the following 2x2 block matrix A^ := - ^ lm,m-2fc 

J-m— 2fc,m 'J 



Then 



haikA 



haf Ai 



[m 



2ky: 



{3.1] 
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Proof. It is enough to consider the nontrivial case /c < y. Let = 
{y{A), Ew{A)), Gw,k = (y{A.k), Ew{Ak)) be the weighted graphs corresponding 
to A, Ak respectively. Note that Gyj is a weighted subgraph of G^i^k induced by 
V{A) = (m) C {2m -2k) = V{Ak). Furthermore, each vertex in V{Ak)\V{A) 
is connected exactly to each vertex in V^(^). The weights of each of these edges 
is 1. These are all edges in G{Ak). A perfect match in G{Ak) correspond to: 

• An m — 2k match between the set of vertices V{Ak)\V{A) and the set 
of vertices a' G Qm-2k,m, viewed as a subset of V{A). 

• A k match between the set of vertices a := (m)\a' C V{B). 

Fix a G Q2k,m- Then the total weight of A;-matchings in Gu,{Ak) using the 
set of vertices a C V{Ak) is given by haf A[q:, a]. The total weight of m — 2k 
matchings using V{Aj^)\V{A) and V{A)\a is {m — 2k)\. Hence the total weight 
of perfect matchings in (^^(A^), which matches the set of vertices a C V{Ak) 
is given by (m — 2k) \ haf A[a, a]. Thus haf Ak = {m — 2k) \ haf^ A. □ 

It is not known if the computation of the number of perfect matching in an 
arbitrary undirected graph on an even number of vertices, or more generally 
the computation of haf A for an arbitrary A G S2n,o(I^+)) has a fpras. The 
probabilistic algorithm outlined in [7] carries over to the computation of haf A, 
however it is an open problem if this algorithm is a fpras. Theorem 13.11 shows 
that the computation of haffc A, for A G Sm_o(I^+)! has the same complexity 
as the computation of haf A, for A G S2n.o(l^+)- 



4 Remarks 

In this section we offer an explanation, using the recent results in why 
perm A is a nicer function than haf A. Let A = [aij] G Sn(M), B = [by] G W^"". 
Forx:= G R" let 

n n ^ 

i=i j=i 

Then permS = q^^'^q^ and haf ^4 = ((f )!)~^ a3;i^"ax ^(^)^ ^ even. 
Assume that B G M"^" has no zero row. Then j9(x) is a positive hyperbolic 
polynomial. (See the definition in [Sj.) Assume that A G S2m,o(R+) is irre- 
ducible. Then g(x), and hence any power q'(x)*,z G N, is positive hyperbohc 
if and only if all the eigenvalues of A, except the Perron-Frobenius eigenvalue, 
are nonpositive. 
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